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Sex determination in skeletal 
remains from the medieval 
Eastern Adriatic coast - 
discriminant function analysis 
of humeri 



Aim To investigate the usefulness of humerus measure- 
ment for sex determination in a sample of medieval skel- 
etons from the Eastern Adriatic Coast. Additional aim was 
to compare the results with contemporary female popula- 
tion. 

Methods Five humerus measurements (maximum length, 
epicondylar width, maximum vertical diameter of the head, 
maximum and minimum diameter of the humerus at mid- 
shaft) for 80 male and 35 female medieval and 19 female 
contemporary humeri were recorded. Only sufficiently pre- 
served skeletons and those with no obvious pathological 
or traumatic changes that could affect the measurements 
were included. For ten samples, analysis of DNA was per- 
formed in order to determine sex using amelogenin. 

Results The initial comparison of men and women indi- 
cated significant differences in all five measures (P< 0.001). 
Discriminant function for sex determination indicated that 
as much as 85% of cases could be properly categorized, 
with better results in men (86%) than women (80%). Fur- 
thermore, the comparison of the medieval and contempo- 
rary women did not show significant difference in any of 
the measured features. Sex results obtained by anthropo- 
logical and DNA analysis matched in all 10 cases. 

Conclusion The results indicate that humerus measure- 
ment in Croatian medieval population may be sufficient to 
determine the sex of the skeleton. Furthermore, it seems 
that secular changes have not substantially affected con- 
temporary population, suggesting that the results of this 
study are transferable to contemporary population as well. 
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Once the skeletal remains are uncovered, anthropologists 
initially aim to reconstruct the biological profile of the per- 
son, which includes sex, age, and height estimation. During 
the reconstruction process, numerous issues may arise, in- 
cluding bone fragmentation and poor preservation of skel- 
etal remains, coupled with the complexity of human skel- 
eton (1,2). Sex determination is one of the first and basic 
steps of assessing the biological profile. Although the anal- 
ysis of DNA is the most reliable method for sex determina- 
tion (3), it is also the most expensive and time consuming 
method, which can also be hindered by local conditions. 
This may especially be true in cases of poor preservation 
of the remains, inhibitors effects, or a small amount of ex- 
tracted DNA from the sample. 

In absence of DNA results, skeletal remains can be used 
to infer subject's sex via two methods, morphological and 
anthropometric. The morphological approach is based 
on the examination of the bones that show the strongest 
sexual dimorphism, principally the skull and the pelvis (4). 
However, this method is not always reliable, especially if 
the skull is fragmented or incomplete. Age can also affect 
the results, especially in elderly women, in which morpho- 
logical characteristics of the skull tend to resemble those 
of men (5). Although morphological methods are very im- 
portant for a preliminary sex assessment, they additionally 
rely on the experience of the examiner and are therefore 
rather subjective and unreliable. 

The second approach is based on anthropometric analysis, 
which relies on the bone measurements. The main analytic 
approach is based on discriminant function analysis, which 
attempts to classify subjects into each of the sexes, by us- 
ing either one or more bones (6). This kind of analysis is a 
very important quantitative method (7) for sex determina- 
tion as it reduces the subjectivity of the examiner (2,4,8). 
So far, only a few such studies have been published using 
Croatian bone samples. These include medieval and con- 
temporary femurs and tibias (8-1 1) and medieval and con- 
temporary mandibles and teeth (6,1 2). Such studies are im- 
portant, since clear differences were observed in different 
populations (13-15), making this a locally-specific feature 
that requires the development of regional standards, ap- 
plicable for local population (1 6). 

Besides already used femurs and tibias, humerus is anoth- 
er long bone from the body that is presumably informa- 
tive for sex determination. This idea was initially derived 
from the empirical investigations of the skeletal remains, 
and further supported by the previously reported sexu- 



al dimorphism of humeri (7,17), even in cases of severe 
bone fragmentation (18). Furthermore, such location-spe- 
cific results may be of interest in modern forensics as well, 
since observed changes in the skeleton marked predomi- 
nantly by the increase in height (1 9), appear to be propor- 
tional, with no indication of sexual dimorphism in ancient 
and modern samples (20). Therefore, the aim of this study 
was to investigate the possibility to determine the sex of 
the subject based on anthropometric analysis of humeri 
measures. 

MATERIALS AND METHODS 

We analyzed a sample of 80 male and 35 female humeri 
from 7 medieval sites from eastern Adriatic coast: Svecurje 
-2estinj (dated in 9/1 1th century) (21), Rizinice (dated 
9/1 0th century), Bijaci Stombrate (dated in 9/1 0th centu- 
ry) (22), Ostrovica Greblje (dated in 9th century) (23), Sopot 
Benkovac (dated in 14/1 5th century), Kamenmost Kaldrma 
(dated in 14/1 5th century) (24), and OtokVuletina rupa - 
Grebcine (dated 17/1 8th century) (25). These graveyards 
were shown by the archeological excavations to be of 
typical Croatian culture, creating a relatively homogenous 
group of subjects involved in this analysis. 

In order to compare these results to the contemporary pop- 
ulation, we additionally used another sample of skeletons, 
dated to the end 19th century and beginning of the 20th 
century. This sample consisted of 1 9 female sceletons from 
the Kozala monastery graveyard (Rijeka, Croatia). The valida- 
tion sample selection was based on the available well pre- 
served skeletons. In both samples, only those with overall 
very good preserved status were included, with the subject 
sex determined using standard anthropological methods 
for sex determination. Sex was determined by examination 
of sex specific pelvis and skull morphological characteris- 
tics. On the pelvis we examined the greater sciatic notch, 
pelvis size and shape, ventral arc, subpubic concavity, and 
medial aspect of ischiopubic ramus, and on the skull: nuchal 
crest, mastoid process, supraorbital margin and ridge, men- 
tal eminence, and ramal flexion (4,26-28). Only skeletons 
meeting the criteria of full confidence in sex determination 
were included in the study. In turn, all of the samples where 
substantial damage was recorded, samples with obvious 
pathological findings, traumas, and deformations that can 
alter the bone structure or affect the measurements were 
excluded from the analysis. A total of 459 medieval adult 
skeletons were analyzed, 1 15 of which met the described 
criteria, while out of 47 contemporary female skeletons, 
1 9 met the described criteria. 
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A total of five measures were made: 

1. maximum length of the humerus - the direct distance 
from the most superior point on the head of the humerus 
to the most inferior point on the trochlea, 

2. epicondylar width of the humerus - the distance of the 
most laterally protruding point on the lateral epicondylefrom 
the corresponding projection of the medial epicondyle, 

3. maximum vertical diameter of the head of the humerus 
- the direct distance between the most superior and infe- 
rior points on the border of the articular surface, 

4. maximum diameter of the humerus at midshaft - the 
maximum diameter of the midshaft measured by turning 
the bone until the maximum diameter is obtained, 

5. minimum diameter of the humerus at midshaft - the 
least diameter of the midshaft (29). 

The measurements were taken using osteometric board 
and sliding caliper. Two measurements were performed in- 
dependently by two authors and later compared. In case 
of discrepancies between the measured values the mea- 
surement was repeated. Measurements were taken from 
the left side whenever possible. 

Analysis of DNA 

In addition, 10 tooth samples from 10 different individu- 
als from the excavation site Otok Vuletina rupa - Grebiine 
were taken for additional confirmation that morphomet- 
ric sex determination is 98% reliable when both pelvis and 
skull are available (4). For these 1 0 samples, sex was deter- 
mined by analysis of the amelogenin gene using protocols 
described previously (30,31). Contamination of DNA was 
prevented using previously described protocol (32). Teeth 
were washed with distilled water and dried. Using circu- 
lar saw type 900 (KaVo Elektrotechnisches Werk, Vertrieb- 
sgesellschaft GmbH, Leutkirch, Germany), the teeth sur- 
face and canals were removed around 2-3 mm in depth. 
The sawdust was washed with distilled water, and the 
teeth were cleaned with 5% commercial bleach. Using liq- 
uid nitrogen the teeth were frozen and then ground into 
powder in a cylinder crushing tool. For each sample, 1 g of 
teeth powder was used and cleaned with EDTA for 5 days. 
DNA was extracted using the 3 ml_ extraction buffer (10 
umol/L Tris, pH 8.0; 100 umol/L NaCI; 50 umol/L EDTA, 
pH 8.0; and 0.5% sodium dodecyl sulfate, SDS) and 



1 50 uL of 20 mg/mL proteinase K (PK) and incubated for 
24-48 hours at 56°C with shaking. The liquefied teeth were 
rinsed with Phenol: Chloroform: Isoamyl Alcohol (25:24:1) 
and centrifuged twice for 10 minutes at 5000 RPM.The top 
layer of liquefied teeth was placed in a new 15-mL tube 
to be mixed with 3 mL n-butanol and centrifuged for 10 
minutes at 5000RPM.The bottom layer was transferred to 
an Amicon tube and repeatedly washed and centrifuged 
with ddH20 for 10 minutes at 2600RPM (at least twice). 
Amplification was performed using AmpFISTR® MiniFiler 
PCR Amplification Kit (Applied Biosystems, Foster City, CA, 
USA) and polymerase chain reaction (PCR) was performed 
in GeneAmp PCR System 9700 (Applied Biosystems). PCR 
products were typed on ABI Prism 310 Genetic Analyzer 
(Applied Biosystems). 

Statistical analysis 

Descriptive statistics encompassed numbers and per- 
centages for categorical variables, coupled with means 
and standard deviations for numerical variables. Inferen- 
tial statistics was based on an independent f test used for 
analysis of sex-related differences in measured variables. 
Discriminant analysis was used to define the existence of 
sex-discriminatory variables, with calculation of percent 
of correctly classified cases as the validation measure. In 
line with previous analyses, we also employed a regression 
analysis, which aimed to identify sex-related differences, in 
line with similar previous studies (6,8,10,1 1,33). All analyses 
were performed using SPSS (ver 18; SPSS Inc, Chicago, IL, 
USA), with the significance level set at P < 0.05. 

RESULTS 

This study was based on 1 1 5 bone samples, 80 male and 35 
female, from 7 medieval locations (Table 1). Most of the an- 
alyzed medieval humeri were from the leftside (80 samples: 
26 women, 54 men), while a smaller portion was from the 
right side (35 samples: 9 women, 24 men). Of the contem- 
porary female samples, 19 were left humeri. The compari- 
son of the left and right humeri did not yield significant dif- 
ference for any of the measured features (data not shown). 

Since all of them were of comparable bone status, we per- 
formed a direct comparison, which suggested that all five 
humerus measures in the medieval period were signifi- 
cantly different between men and women (Table 2). 

The main step of the analysis involved discriminatory 
function, which suggested that the humeri were more ef- 
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fective in sex determination of men than women (Table 
3). In total, 84.8% of the subjects were correctly classified 
using this function. Additionally, we aimed to see if there 
were fixed combinations of the measurement that would 
provide better sex determination functions. When maxi- 
mum and minimum head diameters were used in discrim- 
ination analysis, the results indicated that 75.2% of cases 
were correctly classified. In turn, the use of head width 
and epycondilar width suggested that 82.5% of the sub- 
jects were correctly classified, thus suggesting that all five 
measurements provided the most informative sex deter- 
mination set. 

Furthermore, we aimed to see if it was possible to define 
the range of measurements which provided a-priori cor- 
rect classification for men and women. The results suggest- 
ed that the majority of the measures had a wide range of 
overlap, with the worst result for maximum head diameter, 
which seemed to have the greatest overlap between men 
and women (Table 4). 

Lastly, we aimed to compare the medieval and contem- 
porary samples. This analysis suggested that neither of the 
analyzed measures were significantly different between 
the two subsamples (Table 5). 

Sex was determined by amelogenin analysis for 10 skel- 
etons, 9 male and 1 female. Sex results obtained by DNA 
and morphometric anthropological analysis matched in all 
1 0 cases. 



TABLE 1.The number of humeri from each excavation site 



Site 


Women 


Men 


Total 


Svecurje 


3 


5 


8 


Rizinice 


1 


1 


2 


Bijaci Stombrate 


7 


10 


17 


Ostrovica Greblje 


10 


25 


35 


Sopot Benkovac 


2 


10 


12 


Kamenmost Kaldrma 


7 


10 


17 


Otok Vuletina rupa - Grebcine 


5 


19 


24 


Total 


35 


80 


115 



DISCUSSION 

The results of this study show that humeral measurements 
of the medieval Croatian population may serve as the rea- 
sonably good estimate of sex. As humeri had not before 
been analyzed in the Croatian population for this purpos- 
es, the main aim of this study was to test if the humeral 
measurements were a reliable sex indicator. Determination 
of sex is the first step in determination of biological profile 



TABLE 3. Accuracy for sex determination functions for medi- 
eval period (mean ± standard deviation) 



Measurement; mm 


Men 


Women 


All variables 


58±86.57 


20 ±80.00 


Maximum length 


69 ±90.79 


17± 50.00 


Epicondilar width 


65 ±8784 


19 + 65.52 


Maximum vertical head diameter 


62 ±84.93 


24 ±72.73 


Maximum diameter at midshaft 


75 ±93.75 


11 ±50.00 


Minimum diameter at midshaft 


71 ±88.75 


18± 52.94 


TABLE 4. Overlap ranges of the analyzed measurements for 
men and women 




Sex range 


Measurement; mm 


female 


unclear male 


Maximum length 


<247 


247-385 >385 


Epicondilar width 


<53 


53-65 >65 


Maximum vertical head diameter 




38-46 >46 


Maximum diameter at midshaft 


<20 


20-23 >23 


Minimum diameter at midshaft 


<16 


16-20 >20 



TABLE 5. Comparison of the medieval and contemporary 
women (mean ± standard deviation) 





Medieval 


Contemporary 




Measurement; mm 


women (n = 35) women (n = 1 9) 


P 


Maximum length 


30424 ±17.79 


303.14± 18.62 


0.869 


Epicondilar width 


56.74 ±3.60 


56.1 5 ±4.98 


0.637 


Maximum vertical head 


41 .65 ±2.04 


41.57 + 2.83 


0.984 


diameter 








Maximum diameter at 


21.10+1.45 


20.38 ±2.04 


0.160 


midshaft 








Minimum diameter at 


16.92 ±1.64 


17.23 + 2.13 


0.369 



midshaft 



TABLE 2. Initial comparison of the humeri measurement in men and women from medieval period (mean ± standard deviation) 



Measurement; mm Men(n = 80) Women (n = 35) Index P 



Maximum length 


327.92 ±19.05 


304.24 ±17.79 


107.78 


<0.001 


Epicondilar width 


63.53 ±4.55 


56.74 ±3.60 


111.97 


<0.001 


Maximum vertical head diameter 


46.86 ±3.65 


41 .65 ±2.04 


112.51 


<0.001 


Maximum diameter at midshaft 


23.38 ±1.96 


21.10+ 1.45 


110.81 


<0.001 


Minimum diameter at midshaft 


19.60±2.25 


16.92 ±1.64 


115.84 


<0.001 
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of a person, that is, the first step in individualization, and 
in forensic sciences - identification of an individual. As the 
morphological method of sex determination is subjective 
and relies mostly on the experience of the examiner, an- 
thropometric methods have been developed. These meth- 
ods include discriminant functions for sex determination 
for almost every bone in the human body. But, as reported 
by various authors, these functions tend to be population- 
specific, therefore, the imperative of every region is to de- 
velop its own functions (13-1 6). 

We obtained better classification results for men, and the 
overall pattern suggested that measurements of the en- 
tire humerus provided the best determination possibilities, 
somewhat better than the isolated central parts of humer- 
us or its proximal fragments. These findings are largely in 
line with previous studies (7,1 7,30), which have also point- 
ed out population-specific estimates (1 3-1 6). 

When using one function, the most accurate function 
for women is the maximum head diameter (which clas- 
sifies correctly 72.73% of women), while for men it is the 
maximal diameter at midshaft (which classifies correctly 
93.75%of men). It is interesting that the largest gap be- 
tween men and women is visible in both of these func- 
tions: the most accurate function for male sex determi- 
nation is also the most inaccurate function for female sex 
determination and vice versa. This can be the result of vari- 
ability between sexes, but a difference in a sample size has 
also to be considered. 

Prediction of sex is of higher accuracy in men, which was 
supported by other authors (34). 

The sex difference in the humeral measurements are prob- 
ably due to differential bone remodeling between sexes, in 
men cortical bone develops more during adolescence (35). 

The two measurements with greatest sex difference are the 
maximum length and the maximum diameter at midshaft, 
which was also found by other authors (36). Some authors 
believe that this is common in populations with extreme- 
ly high or extremely low protein consumption (37), while 
other suggest that the circumferential measurements are 
more important for sex determination because of the in- 
fluence of physical activity on bone (38). Some authors 
found that the most effective single dimensions were ver- 
tical head diameter (18,39) and epicondylar breadth (40), 
which indicates the necessity of developing regional 
sex discriminant functions. 



One of the interesting findings is the similarity of the hu- 
meral measurements for medieval and contemporary 
population. This result suggests that any secular changes 
that were happening over time did not affect humeri. An 
important conclusion that can be drawn from this study is 
the possibility to use the results inter-changeably in both 
medieval and contemporary populations. This finding is 
very interesting, since similar studies performed on femurs 
and tibias showed the opposite result (8-1 1 ). Both of these 
conclusions can be joined into a more general result, sug- 
gesting that some parts of the skeleton were seemingly 
more affected than others. An analysis of contemporary 
male humeral samples is necessary to see if the male sam- 
ple will follow the same trend. This is important because 
if the residence is patrilocal, the variability between male 
skeletons will be small inside one cemetery, but consider- 
able between various cemeteries. On the other hand, fe- 
male skeletons will vary more inside one cemetery than 
between cemeteries (41). 

There are several study limitations. First, the overall sample 
size was relatively small, especially the contemporary pop- 
ulation. The analyzed medieval sample encompassed a 
long stretch of time and probably a diverse set of environ- 
mental conditions, which could have affected the final re- 
sult. The results reported here are limited only to well-pre- 
served skeletons, without any substantial skeletal changes 
that could affect the results. 

Ourfindings, showing homogeneity and lack of significant 
differences between medieval and contemporary popula- 
tions are in line with some of the previous studies, which 
suggested homogeneity of the population over time, 
based on the Y chromosome (42) and mitochondrial DNA 
(43). Such findings are at the moment based on limited 
amount of evidence, suggesting that future studies should 
aim to encompass larger periods of time and more pre- 
cisely defined populations in order to establish a solid un- 
derstanding of the anthropometric and genetic changes 
in the population of contemporary Croatia. 
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